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ABSTRACT 

We address the effects of cosmological surface brightness dimming on ob- 
servations of faint galaxies by examining the distribution of "unobscured" star 
formation rate intensities versus redshift. We use the star formation rate in- 
tensity distribution function to assess the ultraviolet luminosity density versus 
redshift, based on our photometry and photometric redshift measurements of 
faint galaxies in the HDF and the HDF-S WFPC2 and NICMOS fields. We find 
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that (1) previous measurements have missed a dominant fraction of the ultravio- 
let luminosity density of the universe at high redshifts by neglecting cosmological 
surface brightness dimming effects, which are important at redshifts larger than 
z ~ 2, (2) the incidence of the highest intensity star forming regions increases 
monotonically with redshift, and (3) the ultraviolet luminosity density plausibly 
increases monotonically with redshift through the highest redshifts observed. By 
measuring the spectrum of the luminosity density versus redshift, we also find 
that (4) previous measurements of the ultraviolet luminosity density at redshifts 
z < 2 must be reduced by a factor 2 to allow for the spectrum of the luminosity 
density between rest-frame wavelengths 1500 and 2800 A. And by comparing with 
observations of high-redshift damped Lya absorption systems detected toward 
background QSOs, we further find that (5) the distribution of star formation 
rate intensities matches the distribution of neutral hydrogen column densities at 
redshifts z ~ 2 through 5, which establishes a quantitative connection between 
high-redshift galaxies and high column density gas and suggests that high-redshift 
damped Lya absorption systems trace lower star formation rate intensity regions 
of the same galaxies detected in star light in the HDF and HDF-S. Because our 
measurements neglect the effects of obscuration by dust, they represent lower 
limits to the total star formation rate density. 

Subject headings: cosmology: observations; galaxies: evolution 

1. INTRODUCTION 

The rest-frame ultraviolet luminosity density of the universe traces the history of cosmic 
star formation because it is produced by hot, massive, young stars. Previous measurements 
have found that the ultraviolet luminosity density increases with redshift to a maximum at 
redshift z ~ 1 to 2 and then decreases (Madau et al. 1996) or remains constant (Madau, 
Pozzetti, & Dickinson 1998; Steidel et al. 1999; Hopkins, Connolly, & Szalay 2000) with 
redshift to higher redshifts. These previous measurements have been interpreted in the 
context of galaxy formation and evolution scenarios as indicating that the stellar content of 
galaxies was formed gradually, over most of cosmic time (e.g. Somerville, Primack, & Faber 
2001). 

But previous measurements have neglected cosmological surface brightness dimming ef- 
fects, which can play a dominant role in setting what is observed at high redshifts (Pascarelle, 
Lanzetta, & Fernandez-Soto 1998). Specifically, surface brightness (per unit frequency in- 
terval) decreases with redshift as (1 + z)~ 3 due to the expansion of the universe. This has an 
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important consequence for observations of distant galaxies: while both intrinsically bright 
and intrinsically faint regions of low-redshift galaxies are accessible to observation, only in- 
trinsically bright regions of high-redshift galaxies are accessible to observation — intrinsically 
faint regions of high-redshift galaxies are of low surface brightness and are simply not de- 
tected against the background noise. Thus all measurements miss some fraction of the light 
of distant galaxies; this fraction is small for low-redshift galaxies but can be dominant for 
high-redshift galaxies. 

Here we address this issue by examining the distribution of "unobscured" star forma- 
tion rate intensities versus redshift — i.e. the star formation rate intensities that are directly 
inferred from the observed rest-frame ultraviolet light that is not obscured by intervening 
dust. We use the star formation rate intensity distribution function to assess the ultraviolet 
luminosity density versus redshift, based on our photometry and photometric redshift mea- 
surements of faint galaxies in the Hubble Deep Field (HDF) and Hubble Deep Field South 
(HDF-S) WFPC2 and NICMOS fields (Lanzetta, Yahil, & Fernandez-Soto 1996; Fernandez- 
Soto, Lanzetta, & Yahil 1999; Yahata et al. 2000; Yahata et al. 2001 in preparation). The 
star formation rate intensity distribution function bears on the issue because it makes ex- 
plicit the effects of cosmological surface brightness dimming on observations of high-redshift 
galaxies, thereby allowing the amount by which previous measurements have underestimated 
the ultraviolet luminosity density at high redshifts to be directly estimated. Because our 
measurements neglect the effects of obscuration by dust, they represent lower limits to the 
total star formation rate density. 

The framework of the unobscured star formation rate intensity and the star formation 
rate intensity distribution function is presented in § 2. The photometry and photometric 
redshift measurements are summarized in § 3, the angular area versus depth relations are 
summarized in § 4, and the Stony Brook faint galaxy redshift survey is summarized in § 5. 
The star formation rate intensity distribution function measurements are described in § 6, 
results are described in § 7, and the connection with the neutral hydrogen column density 
distribution function is described in § 8. The summary and conclusions are presented in § 9. 
Unless otherwise stated, we adopt a standard Friedmann cosmological model of dimensionless 
Hubble constant h wo = H /(100 km s _1 Mpc~ 3 ) and deceleration parameter q = 0.5; this is 
the same cosmological model adopted for previous measurements of the rest-frame ultraviolet 
luminosity density. 
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2. FRAMEWORK 

To address the issue of the effects of cosmological surface brightness dimming, we ex- 
amine the distribution of "unobscured" star formation rate intensities versus redshift. We 
designate the unobscured star formation rate intensity as x and the star formation rate in- 
tensity distribution function as h(x). The star formation rate intensity distribution function 
h(x) bears on the issue because it makes explicit the effects of cosmological surface brightness 
dimming on observations of high-redshift galaxies. 

By "unobscured star formation rate intensity," we mean the star formation rate inten- 
sity that is directly inferred from the rest-frame ultraviolet light that is not obscured by 
intervening dust. Hence unobscured star formation rate intensity (measured, say, in units of 
M Q yr -1 kpc -2 ) is equivalent to rest-frame ultraviolet luminosity per unit area (measured, 
say, in units of erg s _1 Hz -1 kpc~ 2 ), once some fixed scaling between ultraviolet luminosity 
and star formation rate (which is set by the stellar initial mass function) is specified. To 
determine the total star formation rate intensity from the unobscured star formation rate 
intensity would require knowledge of the effects of obscuration by intervening dust, which is 
beyond the scope of the present analysis. Rather, here we address only the effects of cosmo- 
logical surface brightness dimming. We emphasize that because our measurements neglect 
the effects of obscuration by dust, they represent lower limits to the total star formation rate 
density. 

We define the star formation rate intensity distribution function h(x) in such a way 
that h(x)dx is the projected proper area per comoving volume of unobscured star formation 
rate intensity in the interval x to x + dx. The star formation rate intensity distribution 
function h(x) is a fundamental statistical description of the evolving galaxy population, 
similar in spirit to the galaxy luminosity function but cast in terms of star formation rate 
intensity (derived from measurement of surface brightness) rather than luminosity (derived 
from measurement of energy flux). For our purposes, the first moment of h(x) 

oo 

x h(x) dx (1) 

is of particular interest. The first moment of h(x) is the unobscured star formation rate den- 
sity ps, which is equivalent to within a factor of scale to the ultraviolet luminosity density 
p v . Equation (1) makes explicit the effects of cosmological surface brightness dimming on 
measurement of the ultraviolet luminosity density: In equation (1), the limits of integration 
extend from x = through oo. But in practice, a given observation (at a given redshift) 
is sensitive only to some star formation rate intensity threshold x m i n , which is set by the 
surface brightness threshold of the observation (and the redshift). Whether or not a given 
observation (at a given redshift) is suitable for measuring the ultraviolet luminosity density 
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depends on whether or not x min is less than the value of x above which the bulk of the 
ultraviolet luminosity density is emitted. Here we assess the ultraviolet luminosity density 
versus redshift and determine the amount by which previous measurements have underes- 
timated the ultraviolet luminosity density at high redshifts, based on our photometry and 
photometric redshift measurements of faint galaxies in the HDF and HDF-S WFPC2 and 
NICMOS fields. 

3. PHOTOMETRY AND PHOTOMETRIC REDSHIFT MEASUREMENTS 

3.1. Photometry 

Details of our photometry of faint galaxies in the HDF and HDF-S WFPC2 and NIC- 
MOS fields have been described previously by Fernandez- Soto et al. (1999) and Yahata et 
al. (2000). In this section, we summarize important aspects of our photometry techniques. 

We processed all available space- and ground-based optical- and infrared-wavelength 
images of the HDF and HDF-S WFPC2 and NICMOS fields. These images include discre- 
tionary, GO, and GTO observations obtained with the Hubble Space Telescope (HST) and 
images obtained with the Kitt Peak National Observatory (KPNO) 4 m telescope and the 
European Southern Observatory (ESO) New Technology Telescope (NTT) and Very Large 
Telescope UT-1 (VLT UT-1). We registered the images (to an accuracy of typically 0.25 
pixel for space-based images and 0.5 pixel for ground-based images) by appropriate trans- 
lation, rotation, and scaling using our own algorithms and software, and we identified and 
masked cosmic ray events using our own algorithms and software. Details of the observations 
are summarized in Table 1. 

We detected objects in the images at multiple bandpasses using the SExtractor pro- 
gram (Bertin & Arnouts 1996), starting at shorter wavelengths and working toward longer 
wavelengths. First, we detected objects in the images of a fiducial bandpass (typically the 
F814W bandpass). Then, we masked regions around detected objects and detected objects 
in the unmasked regions of images of a longer-wavelength bandpass (typically the F160W 
or H bandpass). Finally, we repeated these steps for remaining images of longer- wavelength 
bandpasses (typically the F222M or K bandpass). In this way, the shorter- wavelength band- 
passes (which are generally of higher sensitivity and resolution but which are unsuitable 
for detecting high-redshift galaxies) were used to detect low- and moderate-redshift galaxies 
and the longer- wavelength bandpasses (which are generally of lower sensitivity and resolution 
but which are suitable for detecting high-redshift galaxies) were used to detect high-redshift 
galaxies. Details of the object detection are summarized in Table 2. 
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We performed photometry using our quasi-optimal photometry technique, which fits 
model spatial profiles of the detected objects to the space- and ground-based images (Yahata 
et al. 2000). We determined model spatial profiles by using non-negative least squares 
image reconstruction (Puetter & Yahil 1999) of one or more of the space-based images. 
The quasi-optimal photometry technique offers three important advantages in comparison 
with conventional methods: (1) it provides higher signal-to-noise ratio measurements, (2) 
it accounts for point-spread function variations between the images, and (3) it accounts for 
uncertainty correlations between nearby, overlapping neighbors. 



3.2. Photometric Redshift Measurements 

Details of our photometric redshift measurements of faint galaxies in the HDF and 
HDF-S WFPC2 and NICMOS fields have been described previously by Lanzetta et al. 
(1996), Fernandez-Soto et al. (1999), and Yahata et al. (2000). In this section, we summarize 
important aspects of our photometric redshift measurement techniques. 

We measured photometric redshifts (and spectral types) of galaxies using our redshift 
likelihood technique (Lanzetta et al. 1996; Fernandez-Soto et al. 1999) with a sequence of 
six spectrophotometric templates, including templates of E/S0, Sbc, Scd, and Irr galaxies 
and low- and high-extinction starburst galaxies (which we designate as SB1 and SB2), and 
incorporating the effects of intrinsic and intervening absorption by neutral hydrogen. 

We assessed the accuracy and reliability of the photometric redshift measurements by 
comparing photometric redshift measurements with spectroscopic redshift measurements of 
galaxies identified in the HDF and HDF-S. Results indicate that at all redshifts z < 6 that 
have yet been examined, the photometric redshift measurements are characterized by an RMS 
relative dispersion with respect to the spectroscopic redshift measurements of Az/(1 + z) < 
0.065 and that there are no known examples of photometric redshift measurements that 
differ from spectroscopic redshift measurements by more than a few times the RMS relative 
dispersion (Yahata et al. 2000; Fernandez-Soto et al. 2001). These results set an upper 
limit to the systematic uncertainty of the photometric redshift technique, to which must be 
added, of course, the uncertainty due to effects of photometric error to determine the total 
uncertainty of a particular measurement. 
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4. ANGULAR AREA VERSUS DEPTH RELATION 

The angular area versus depth relations set the angular area accessible to the observa- 
tions as a function of "depth," where "depth" might in general be any one of a number of 
observed or intrinsic properties of galaxies (e.g., energy flux in a given observed-frame band- 
pass, luminosity at a given rest-frame wavelength, or surface brightness at a given rest-frame 
wavelength). For the purpose of the present analysis, we take "depth" to be unobscured star 
formation rate intensity x. The angular area versus depth relation must account for several 
important effects: First, the sensitivities of the individual images (especially the NICMOS 
images) vary with position. Second, objects are detected in one or more observed-frame 
bandpasses, whereas results are determined with respect to a single observed-frame band- 
pass or rest-frame wavelength. Finally, object detection is performed at multiple bandpasses, 
which introduces complicated conditional selection criteria into the analysis. 

We determined the angular area versus depth relation as a function of redshift z and 
unobscured star formation rate intensity x using surface brightness sensitivity maps of the 
individual HDF and HDF-S WFPC2 and NICMOS images. First, we determined the surface 
brightness sensitivity maps by scaling the sky variance maps in such a way as to match 
the faint envelopes of surface brightnesses actually included into the segmentation maps 
of the galaxies. (This assumes that the sky variance maps trace the shape but not the 
normalization of the surface brightness sensitivity maps.) Then, we calculated the surface 
brightness expected for given values of z and x, using the K correction of an Irr galaxy. (No 
cosmological model is required, because the dependence of surface brightness on redshift is 
independent of cosmological model.) Then, we summed the angular areas of the pixels that 
are sensitive enough to detect the surface brightness expected for the given values of z and 
x, which we identified by examining the surface brightness sensitivity maps of the images 
used for object detection. Finally, we repeated these steps for a range of values of z and x. 

The angular area versus depth relation fl(z, function of redshift z and star 

formation rate intensity x is shown in Figure 1. In Figure 1, white regions represent Q(z, x) = 
0, black regions represent Q(z,x) = 11.7 arcmin 2 (i.e. the total angular area spanned by the 
observations), and grey regions represent intermediate values of Q(z,x). 

5. STONY BROOK FAINT GALAXY REDSHIFT SURVEY 

Our catalogs of photometry and photometric redshift measurements include nearly 5000 
faint galaxies, of which nearly 1000 yield redshift measurements z > 2 and more than 50 yield 
redshift measurements z > 5 (ranging up to and beyond z = 10). The catalogs of photometry 
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and photometric redshift measurements and the angular area versus depth relations together 
constitute a galaxy redshift survey to the faintest limits and highest redshifts yet accessible, 
which we designate as the Stony Brook faint galaxy redshift survey. 

Properties of the survey are summarized in Figure 2, which shows the redshift measure- 
ment distributions of all galaxies identified in the HDF and HDF-S WFPC2 and NICMOS 
fields. The redshift measurement distributions of galaxies identified in the HDF and HDF-S 
WFPC2 field are characterized by broad peaks at redshift z ~ 1 and tails extending to 
redshifts z > 5. (The distributions are statistically different from one another, with the 
HDF-S WFPC2 field exhibiting a statistically significant excess of galaxies of redshift z > 2 
compared with the HDF, and both distributions exhibit statistically significant large-scale 
fluctuations. Because the HDF-S fields were chosen due to their proximity to a known QSO 
of redshift z = 2.2, the HDF-S field may be biased toward galaxies of redshift z « 2.) 
The redshift measurement distribution of galaxies identified in the HDF-S NICMOS field is 
characterized by a broad peak at redshift z « 1 and a tail extending to redshifts z > 10. 

We note several points about the nature of the galaxies of redshift measurement z > 
6 identified in our survey. First, our photometric redshift measurements have been di- 
rectly compared with spectroscopic redshift measurements at redshifts z < 6 but have not 
been directly compared with spectroscopic redshift measurements at redshifts z ~> 6. (See 
Fernandez-Soto et al. 2001 for the latest description of this comparison and for an assessment 
of the accuracy and reliability of the photometric and spectroscopic redshift measurement 
techniques.) Second, the redshift likelihood functions of most (but not all) of the galaxies 
of redshift measurement z > 6 are sufficiently broad or multiply-peaked that lower-redshift 
solutions cannot be excluded. Hence a "best-fit" redshift measurement z > 6 does not 
by itself indicate a redshift z > 6 — in every case the redshift likelihood function must be 
consulted to identify the range of allowed redshifts. And third, our bootstrap error analysis 
(described below) explicity accounts for cases of broad or multiply-peaked redshift likelihood 
functions. The emphasis of the present analysis is on redshifts z < 6, and the galaxies of 
redshift measurement z > 6 do not bear significantly on our primary conclusions. 



6. STAR FORMATION RATE INTENSITY DISTRIBUTION FUNCTION 

MEASUREMENTS 

The objective of the analysis is to determine the star formation rate intensity distribution 
function by considering all pixels in all galaxies in all fields on an individual pixel-by-pixel 
basis. In this section, we describe our measurements of the star formation rate intensity 
distribution function. 
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First, we determined the rest-frame 1500 A luminosity of each pixel, using the measured 
energy flux of the pixel, the redshift of the pixel (which is set by the redshift of the host 
galaxy), the empirical K correction of the pixel (which is set by the spectrophotometric 
template of the host galaxy), and the cosmological model. Then, we determined the star 
formation rate of each pixel, using a Salpeter stellar initial mass function to scale rest- 
frame 1500 A luminosity to star formation rate. Then, we determined the proper area of 
each pixel, using the redshift of the pixel, the angular plate scale of the image, and the 
cosmological model. Then, we divided the star formation rate by the proper area of each 
pixel to determine the star formation rate intensity x of the pixel. The star formation rate 
intensity x is measured, say, in units of M® yr^ 1 kpc~ 2 . 

Next, we summed the proper areas of all pixels within given star formation rate intensity 
and redshift intervals. Then, we determined the comoving volumes that correspond to the 
redshift intervals, using the angular area versus depth relation and the cosmological model. 
Then, we divided the summed proper areas by the appropriate star formation rate intensity 
intervals and comoving volumes to determine the star formation rate intensity distribution 
function h(x) at values of the star formation rate intensity x down to the noise thresholds 
of the images. Finally, we excluded measurements of h(x) at values of x less than the values 
that could have been detected at the 5<r level over 95% of the angular area of the survey. 
The segmentation maps that define the isophotal apertures of the galaxies contain, of course, 
pixels of insignificant or even negative measured energy flux (because the segmentation maps 
are defined with respect to smoothed versions of the images.) Thus low- intensity pixels near 
the edges of the segmentation maps are unsuitable for incorporation into the measurement 
of h(x), since roughly half of these pixels have negative measured energy flux and so cannot 
bear on the analysis. Considering only high-intensity values of x eliminates the possibility 
that low-intensity pixels near the edges of the segmentation maps bias the measurement 
of h(x). We verified that the measurement of h(x) at our adopted star formation rate 
intensity thresholds is not affected by pixel-to-pixel noise by redetermining h(x) from within 
the segmentation maps but with the actual images replaced by Gaussian noise (according 
to the actual noise properties of the images). The star formation rate intensity distribution 
function h(x) is measured, say, in units of kpc 2 (M yr" 1 kpc~ 2 )~ x Mpc~ 3 . 

We determined measurement uncertainties of the star formation rate intensity distribu- 
tion function h(x) by applying a variation of the "bootstrap" resampling technique described 
previously by Lanzetta, Fernandez- Soto, & Yahil (1998). The bootstrap resampling tech- 
nique explicitly accounts for effects of sampling uncertainty, photometric uncertainty, and 
variance with respect to the spectrophotometric templates. First, we resampled objects from 
the catalogs (allowing the possibility of duplication), added random noise to the photometry 
(according to the actual noise properties of the images) and redetermined the photometric 
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redshift measurements, and added random noise to the photometric redshift measurements 
(according to the actual noise properties of the photometric redshift technique, as determined 
from the comparison of photometric and spectroscopic redshift measurements). Then, we 
determined the star formation rate intensity distribution function h(x) from the resampled, 
perturbed catalogs of photometry and photometric redshifts exactly as from the actual cat- 
alogs of photometry and photometric redshifts. Finally, we repeated these steps a large 
number of times in order to trace out the distributions of measurement uncertainties. 



7. RESULTS 

Results of our analysis are shown in Figure 3, which plots the star formation rate 
intensity distribution function h(x) versus the star formation rate intensity x as a function 
of redshift z. Several points are evident on the basis of Figure 3: 

First, the star formation rate intensity threshold of the observations is a steep function 
of redshift, ranging from x ~ 5 x 10~ 4 M yr" 1 kpc" 2 at z — 0.5, to ~ 0.3 M yr" 1 kpc" 2 
at z — 5, to 1 M G yr" 1 kpc" 2 at z — 10. 

Second, at the redshifts z < 1 at which both low- and high-intensity star formation 
rate intensities can be measured, \ogh(x) versus logo; is characterized by a relatively shallow 
dependence at x < 10~ 2 M yr" 1 kpc" 2 and a relatively steep dependence at x > 10" 2 M 
yr" 1 kpc" 2 . The shallow dependence is sufficiently shallow and the steep dependence is 
sufficiently steep that, from equation (1), the bulk of the ultraviolet luminosity density is 
emitted where the two regions join, i.e. at x ps 10" 2 M Q yr" 1 kpc" 2 . This implies that to 
measure the bulk of the ultraviolet luminosity density requires sensitivity to star formation 
rate intensities less than the intensity x ~ 10~ 2 M yr" 1 kpc" 2 at which h(x) steepens. From 
Figure 3, it is clear that the observations of the Hubble Deep Fields — which are the most 
sensitive observations of faint galaxies ever obtained — can detect the intensities at which 
h(x) steepens only to redshift z ~ 1.5 or at most z ~ 2; at higher redshifts, the observations 
are not sensitive to the intensities at which the bulk of the ultraviolet luminosity density are 
emitted. 

Third, the high-intensity end of h(x) — which is directly measured at all redshifts through 
z ~ 10 — evolves more or less monotonically toward higher values at higher redshifts. For 
example, at x ~ 3 M yr _1 kpc~ 2 , h(x) evolves from h(x) w 10~ 5 foioo kpc 2 (M Q yr _1 
kpc~ 2 ) -1 Mpc~ 3 at z w 1.5 through ps 10~ 4 h 100 kpc 2 (M yr" 1 kpc" 2 )" 1 Mpc -3 at z ps 2.5 
through ~ 10" 3 h W o kpc 2 (M yr -1 kpc" 2 )" 1 Mpc" 3 at z ps 4. There are no local (say 
z < 1) analogs of the very high intensity (say x > 1 M Q yr" 1 kpc" 2 ) regions that are seen 
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at higher redshifts, at least over the volumes probed by the Hubble Deep Fields. 

To quantitatively establish the effect of the star formation rate intensity threshold of the 
observations on the determination of the star formation rate density, we fitted h(x) versus 
x by a broken double power law model 

h/Vt _ / A(x/x )~ ai (x < x ) , , 

h[X) -\A(x/x o y^ (x>x ), [Z) 

which is described by slopes a.\ and a 2 at low and high intensities, characteristic intensity 
Xo, and normalization A. First, we established fiducial parameters by fitting the model to 
the observations over the redshift interval z = 0.5 to 1. (These redshifts are low enough 
that the observations are sensitive to low values of x but not so low that redshift measure- 
ment uncertainties can significantly affect measurements of luminosity or star formation rate 
intensity.) Results indicate a low-intensity slope 



a high-intensity slope 



«i = -i-44±8:S, (3) 

a 2 = -3.05t°;i, (4) 



and a characteristic intensity 

log x = -1.96 ±0.13 M yr kpc 2 . (5) 

The fiducial model is shown in Figure 3a. Next, we fitted simple variants of the fiducial 
model over other redshift intervals, fixing a± and a 2 and either (1) scaling vertically (i.e. 
fixing Xq and varying A), (2) scaling horizontally (i.e. fixing A and varying x ), and (3) 
scaling the break intensity [i.e. varying Xq and A in such a way that h(x) at low intensities 
does not change]. The variants of the fiducial model are shown in Figure 3b. Finally, we 
determined the star formation rate density versus redshift by integrating over the models, 
to intensities as low as the lowest star formation rate intensity x = 6 x 10~ 4 M© yr -1 kpc~ 2 
observed in the local universe (Kennicutt 1998). 

Results are shown in Figure 4a, which plots the unobscured star formation rate densities 
versus redshift obtained by integrating over the models. For comparison, Figure 4a also plots 
the unobscured star formation rate density versus redshift obtained by Madau et al. (1996, 
1998) by summing over the observed galaxy energy fluxes. From Figure 4a, it is evident that 
the ultraviolet luminosity density obtained by integrating over the models increases more 
or less monotonically with redshift, for each of the three models. Over the redshift interval 
z = through 10, the increase ranges from a factor 10 (for scaling the break intensity of 
the fiducial model) through a factor 250 (for scaling the fiducial model vertically). The 
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ultraviolet luminosity density obtained by Madau et al. (1996, 1998) by summing over the 
observed galaxy energy fluxes significantly underestimates the unobscured star formation rate 
density obtained by integrating over the models at redshifts larger than z ~ 2. The increase 
of the ultraviolet luminosity density with redshift is a direct consequence of the monotonic 
increase of the incidence of the highest intensity star forming regions with redshift and must 
hold unless the actual functional form of h(x) at redshifts z > 2 differs dramatically from 
the functional forms of the three simple variants of the fiducial model. We conclude that the 
unobscured star formation rate density — or equivalently the ultraviolet luminosity density— 
plausibly increases monotonically with redshift through the highest redshifts observed. 

But it also appears from Figure 4a that our measurements and the measurements of 
Madau et al. (1996, 1998) of the unobscured star formation rate density are discordant at 
redshifts z = to 2. This is a consequence not of cosmological surface brightness dimming 
effects but rather of a mismatch of the rest-frame wavelengths from which the star formation 
rate densities were measured. We measured star formation rate densities by scaling from 
measured luminosity densities at rest-frame wavelength 1500 A at all redshifts. In contrast, 
Madau et al. (1996, 1998) measured star formation rate densities by scaling from measured 
luminosity densities at rest-frame wavelength 2800 A at redshifts z < 2 (based on results of 
Lilly et al. 1996 and Connolly et al. 1997) and at rest-frame wavelength 1500 A at redshifts 
z > 2 (based on results from the Hubble Deep Field), adopting a flat f v spectrum to scale 
luminosity density to star formation rate density. Thus our measurements and the measure- 
ments of Madau et al. (1996, 1998) should be consistent only if the rest-frame ultraviolet 
luminosity density of the universe exhibits a flat f u spectrum at redshifts z = to 2. 

To test this possibility, we measured the spectrum of the luminosity density versus 
redshift, based on galaxies observed in the HDF and HDF-S using our catalog of photometry 
and photometric redshifts. Results are shown in Figure 5. The ratio of the luminosity density 
at rest-frame wavelength 2800 A to the luminosity density at rest-frame wavelength 1500 A 
is ~ 2 at redshifts z < 2. In other words, the flat /„ spectrum adopted by Madau et al. 
(1996, 1998) is inconsistent with the measured spectrum of the luminosity density at redshifts 
z < 2. Consequently, the measurements of Madau et al. (1996, 1998) at redshifts z < 2 must 
be reduced by a factor R3 2 to bring them into consistency with measurements of the star 
formation rate density based on rest-frame 1500 A luminosity density. The results corrected 
in this way are shown in Figure 4b. It appears from Figure 4b that our measurements and 
the corrected measurements of Madau et al. (1996, 1998) of the unobscured star formation 
rate density are consistent at redshifts z < 2. 

Adoption of the currently favored non-zero cosmological constant cosmological model 
of vaccuum energy density fl\ = 0.7 and matter density VLm = 0.3 does not qualitatively 
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change our conclusions. Under this cosmological model, the values of the unobscured star 
formation densities are reduced, by amounts that range from m 20% for the lowest-redshift 
point to i=s 50% for the highest-redshift point. 

8. CONNECTION WITH THE NEUTRAL HYDROGEN COLUMN 
DENSITY DISTRIBUTION FUNCTION 

The star formation rate intensity distribution function h(x) is exactly analogous to the 
column density distribution function f(N) (as a function of neutral hydrogen column den- 
sity N) measured from damped Lja absorption systems detetected toward high-redshift 
background QSOs in that h(x)dx and f(N)dN both measure incidence per length of cosmo- 
logically distributed material. The connection between h(x) and f(N) can be made explict 
by adopting a "Kennicutt" (1998) relation between star formation rate intensity x and neu- 
tral hydrogen column density N 



(Kennicutt 1998 multiplied observed energy fluxes by a factor 2.8 to correct for effects of 
obscuration by intervening dust. Here we divide the energy fluxes of Kennicutt 1998 by 
this factor 2.8 to uncorrect for effects of obscuration by intervening dust, i.e. to establish 
a relation between unobscured star formation rate intensity x and gas column density N. 
We consider neutral hydrogen column density to represent gas column density because the 
molecular content of high column density QSO absorption systems is measured to be low.) 
Then h(x) and f(N) are related as 



where dl is a co moving length element and dX is an absorption distance element. We 
applied equations (6) and (7) to express measurements of f(N) at redshifts z ~ 1.6 through 
5 (Lanzetta et al. 1991; Storrie-Lombardi, Irwin, & McMahon 1996) in terms of h(x). 

Results are shown in Figure 3b. The amplitude of the measurements of f(N) at redshifts 
z ~ 1.6 through 5 expressed in terms of h(x) match the amplitude of the measurements 
of h(x). Further, application of equation (5) to a power-law form of the column density 
distribution function f(N) oc N" 13 of slope (3 = 1.48 ± 0.30 (Storrie-Lombardi, Irwin, & 
McMahon 1996) predicts a power-law form of the star formation rate intensity distribution 
function h(x) oc x~ ai of slope a.\ = 1.34 ± 0.22, which is in excellent agreement with the 
slope «! = 1.44t°;og measured at redshifts z = 0.5 to I. We take the remarkable agreement 





h(x) dx dl = f(N) dN dX, 



(7) 
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between f(N) expressed in terms of h(x) and h(x) — in terms of both amplitude and slope- 
to (1) suggest that the actual functional form of h(x) at redshifts z > 2 does not differ 
dramatically from the function forms of the three simple variants of the fiducial model (and 
that the third variant is slightly favored over the other two variants) and (2) establish a 
quantitative connection between high-redshift galaxies and high column density gas and 
suggest that high-redshift damped Lja absorption systems trace lower star formation rate 
intensity regions of the same galaxies detected in star light in the HDF and HDF-S. 

Of course, f(N) and h(x) are measured on very different physical scales. Specifically, 
f(N) is measured on a physical scale of ~ 1 pc (which is set by the size of the continuum 
emitting regions of QSOs), whereas h(x) is measured on a physical scale of ~ 200 pc (which 
is set by the angular resolution of the HDF and HDF-S images). The remarkable agreement 
between f(N) expressed in terms of h(x) and h(x) nevertheless establishes that faint galaxies 
in the HDF and HDF-S obey a statistical Kennicutt (1998) relation that is consistent with 
the Kennicutt (1998) relation of nearby galaxies. In other words, the relationship between 
column density measured on a scale of ~ 1 pc and star formation rate intensity measured 
on a scale of ~ 200 pc for faint galaxies in the HDF and HDF-S is on average consistent 
with the relationship between column density and star formation rate intensity measured on 
galactic scales for nearby galaxies. 



9. SUMMARY AND CONCLUSIONS 

We address the effects of cosmological surface brightness dimming on observations of 
faint galaxies by examining the distribution of "unobscured" star formation rate intensities 
versus redshift, which we use to assess the ultraviolet luminosity density versus redshift, 
using our photometry and photometric redshift measurements of faint galaxies in the HDF 
and the HDF-S WFPC2 and NICMOS fields. We find the following results: 

1. Previous measurements have missed a dominant fraction of the ultraviolet luminosity 
density of the universe at high redshifts by neglecting cosmological surface brightness 
dimming effects, which are important at redshifts larger than z ~ 2. 

2. The incidence of the highest intensity star forming regions increases monotonically 
with redshift. 

3. The ultraviolet luminosity density plausibly increases monotonically with redshift through 
the highest redshifts observed. 

4. Previous measurements of the ultraviolet luminosity density at redshifts z < 2 must 
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be reduced by a factor R3 2 to allow for the spectrum of the luminosity density between 
rest-frame wavelengths 1500 and 2800 A. 

5. The distribution of star formation rate intensities matches the distribution of neutral 
hydrogen column densities at redshifts z ~ 2 through 5, which establishes a quan- 
titative connection between high-redshift galaxies and high column density gas and 
suggests that high-redshift damped Lya absorption systems trace lower star forma- 
tion rate intensity regions of the same galaxies detected in star light in the HDF and 
HDF-S. 

Because our measurements neglect the effects of obscuration by dust, they represent lower 
limits to the total star formation rate density. Our analysis suggest that star formation in the 
very early universe may have occurred at a much higher rate than is generally believed and 
that cosmological surface brightness dimming effects cannot be ignored when interpreting 
statistical properties of the high-redshift galaxy population. 

This research was supported by NASA grant NAGW-4422 and NSF grant AST-9624216. 
AFS was supported by a European Union Marie Curie Fellowship. The authors acknowledge 
helpful comments from an anonymous referee. 
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TABLE 1 

DETAILS OF THE OBSERVATIONS 

Field Bandpass Reference 



HDF F300W 1 

F450W 1 

F606W 1 

F814W 1,2 

F110W 3,4 

F160W 3,4 

J 5 

H 5 

K 5 

HDF-S WFPC2 F300W 6 

F450W 6 

F606W 6 

F814W 6 

U 7 

B 7 

V 7 
R 7 
I 7 
J 7 
H 7 
K 7 

HDF-S NICMOS F110W 8 

F160W 8 

F222M 8 

U 9 

B 9 

V 9 
R 9 
/ 9 



REFERENCES— (1) Williams et al. (1996); 
(2) Gilliland, Nugent, & Phillips (1999); (3) 
Thompson et al. (1999); (4) Dickinson et al. 
(2000); (5) Dickinson et al. (2001), in prepara- 
tion; (6) Casertano et al. (2000); (7) da Costa 
et al. (2001); (8) Williams et al. (2000); (9) 
ESO VLT-UT1 Science Verification (1998). 
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TABLE 2 

DETAILS OF THE OBJECT DETECTION 





Field 


Object Detection Bandpasses 


HDF 




F814W, F160W, K 


HDF 


-S WFPC2 


F814W, H, K 


HDF- 


-S NICMOS 


. . F160W, F222M 
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Fig. 1. — Angular area Q accessible to the observations as a function of logarithm of star 
formation rate intensity x and redshift z. Star formation rate intensity x is measured in 
units of M Q yr— 1 kpc~2. White represents Q = 0, black represents Q = 11.7 arcmin 2 , and 
grey represents intermediate values of Q. 

Fig. 2. — Redshift measurement distributions of all galaxies identified in the (a) HDF and 
HDF-S (b) WFPC2 and (c) NICMOS fields. 

Fig. 3. — Star formation rate intensity distribution function h(x) versus star formation rate 
intensity x at redshifts ranging from z = to 0.5 through z = 6 to 10. Vertical error bars 
show la uncertainties (determined by a bootstrap resampling technique), and horizontal 
error bars show bin sizes. (For many points, the vertical error bars are smaller than the 
points.) (a) Black circles show h(x) determined from observations of faint galaxies. Green 
line segments show the best-fit broken double power law model fitted to the observations 
at redshifts z = 0.5 to 1. (b) Black circles show h(x) determined from observations of 
faint galaxies. Magenta open circles show h(x) determined by expressing measurements of 
f(N) in terms of h(x). Line segments show variants of the best-fit broken double power law 
model obtained by scaling vertically (red), scaling horizontally (blue), and scaling the break 
intensity (green). 

Fig. 4. — (a) Unobscured star formation rate densities ps versus redshift z. Points plot 
unobscured star formation rate densities versus redshift obtained by integrating over the 
models, fixing a\ and «2 and either (1) scaling vertically (red squares), (2) scaling horizon- 
tally (blue triangles), and (3) scaling the break intensity (green circles) and the unobscured 
star formation rate density versus redshift obtained by Madau et al. (1996, 1998) by sum- 
ming over the observed galaxy energy fluxes (black open circles). Vertical error bars show 
la uncertainties, which for the red squares, blue triangles, and green circles are determined 
by a bootstrap resampling technique, (b) Same as (a), except that the measurements of 
Madau et al. (1996, 1998) at redshifts z < 2 are reduced by a factor 2 to bring them into 
consistency with measurements of the star formation rate density based on rest-frame 1500 
A luminosity density. 

Fig. 5. — Ratio of the luminosity density at rest-frame wavelength 2800 A to the luminosity 
density at rest-frame wavelength 1500 A at redshifts z < 10, based on galaxies observed in 
the HDF and HDF-S using our catalog of photometry and photometric redshifts. 
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